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Semiconductor  lasers  (SLs)  have  a  wide  range  of  applications  but  a  weak 
optical  feedback  is  enough  to  destabilize  its  normal  output.  These  instabili¬ 
ties  are  generating  higher  intensity  or  frequency  noise  which  are  undesirable 
for  certain  applications.  In  close  collaboration  with  the  Nonlinear  Optics 
group  at  Kirtland  AFB,  Dr.  Erneux  has  proposed  to  apply  new  asymp¬ 
totic  techniques  for  solving  the  laser  equations.  These  techniques  allow  us  to 
determine  analytical  expressions  for  all  bifurcation  points. 

The  simplest  problem  that  models  the  interactions  between  a  SL  and  the 
outside  world  is  the  laser  subject  to  optical  injection.  It  is  formulated  by 
three  equations  for  the  amplitude  of  the  electrical  field,  jEj  —  1  -t-  a,  the 
phase  difference  between  master  and  slave  electrical  fields,  i/i,  and  the  carrier 
density  above  threshold,  Z 

^  =  Z{1  + a)  ■¥r)cos{i)),  •  (1) 

^  =  (2) 


=  -Z(l  +  2P(1  +  a)2)  -  P{2a  -h  a^). 


(3) 


In  these  equations,  time  is  5  =  t/r^  where  Tp  is  the  photon  lifetime.  Cl  denotes 
the  detuning  between  the  frequency  of  the  master  laser  and  the  free-running 
slave  laser.  P  is  the  pumping  current  above  threshold.  77  >  0  is  the  injection 
rate.  T  =  Ta/Tp  is  defined  as  the  ratio  of  the  carrier  and  photon  lifetimes. 
b  is  the  linewidth  enhancement  factor  b.  Typical  values  of  P,  Cl,rf,  T  and  b 
are  P  ~  0.1  —  1,  fl  0.01,  77  ~  0.01,  T  10^  and  &  ~  4  —  6.  The  small 
value  of  77  implies  that  the  amplitude  of  the  electrical  field  is  dose  to  its  value 
at  77  =  0  (i.e.,  \E]  =  1)  and  motivates  the  decomposition  l£?|  =  1  -I-  o.  The 


basic  time  of  the  laser  oscillations  is  the  laser  relaxation  oscillations  frequency 
U  ^  y/2P/T. 

Recent  experiments  have  suggested  that  there  exists  a  domain  of  bistabil¬ 


ity  at  large  negative  values  of  the  detuning.  This  bistability  corresponds  to 
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the  coexistence  of  stable  pulsating  intensity  refines  and  stable  steady  states. 
The  main  objective  of  our  analysis  was  to  understand  analjrtically  the  mecha¬ 
nism  leading  to  this  coexistence.  To  this  end,  we  have  considered  the  limit  of 
large  detunings  (specifically,  the  limit  |A|  =  ]n/(i;|  =  0{b),  A  =  r}/u)  =  0{b) 
as  6  — » CX3  of  Eqs.  (l)-(3)). 

At  a  jfixed  negative  A  and  changing  A  from  zero,  we  find  that  this  domain 
of  bistability  is  bounded  below  by  a  steady  state  limit  point  and  above  by 
a  bifurcation  point  to  quasiperiodic  oscillations  located  at  A  =  Agp  defined 


by  [1] 


Agr  =  (A  <  0). 


where  ^  «  1  is  defined  as  the  damping  rate  of  the  laser 

relaxation  oscillations.  Aqp  is  a  bifurcation  point  from  the  multi-wave  mixing 
regime  to  quasiperiodic  oscillations.  In  the  second  part  of  our  activities  at 
Kirtland,  we  investigate  how  Agp{A)  behaves  as  A  is  progressively  increased. 
We  analyze  Eqs.  (l)-(3)  by  a  new  perturbation  analysis  valid  for  arbitrary 
A  =  0(1)  but  small  A.  Away  from  the  resonance  points  A  =  0,  ±1,  ±2, ..,  the 
basic  multi-wave  mixing  regime  (frequency  <T  “=  jA])  exhibits  a  bifurcation  to 
quasiperiodic  oscillations  (frequencies  <Ti  ==  |A|  and  ca  1)  which  is  located 
at  _ 

Agp  s  J*|a(1-A2)  (A{1  -  A^)  >  o)  .  (5) 


Note  that  (5)  is  matching  the  expression  (4)  for  large  negative  A  and  thus 
(5)  is  the  natural  continuation  of  (4)  as  |A|  becomes  0(1).  Near  the  res¬ 
onant  point  A  =  0,  ±1,  ±2,  ...  special  analyses  are  needed  because  the 
quasiperiodic  oscillatioons  may  lock  into  multi-periodic  oscillations.  Figure  1 
represents  the  stability  diagram  of  the  steady  states  in  the  A  vs  A  diagram 
(the  values  of  the  parameters  are  P  =  0.5,  5  =  5  and  T  -=  1000).  The  locking 
domain  corresponds  to  the  domain  A  >  Ax,p  (Ax,p  is  defined  as  the  limit 
point  of  the  steady  states).  The  steady  states  are  however  unstable  in  the 
re^on  bounded  by  the  curve  A  ==  A^  (A^  is  defined  as  a  stable  Hopf  bifurca¬ 
tion  point).  The  domain  below  the  two  straight  lines  A  =  Alp  correspond  to 
the  multi-wave  mixing  regime.  Both  the  multi-wave  mixing  solution  and  the 
Hopf  bifurcation  solution  lead  to  pulsating  intensity  oscillations  but  are  char¬ 
acterized  by  an  unboimded  and  a  bounded  phase  tp,  respectively.  A  =  Aqp 
corresponds  to  a  bifurcation  point  from  liie  multi-wave  mixing  regime  to 
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quasiperiodic  oscillations  and  is  defined  by  (5). 

We  have  considered  both  a  case  of  negative  A  and  a  case  of  positive 
A.  According  to  our  analysis,  a  bifurcation  to  quasiperiodic  oscillations 
should  be  observed  for  A  slightly  less  than  --1  for  A  not  too  large.  Figure 
2  shows  the  numerical  bifurcation  diagram  of  the  periodic  and  quasiperi¬ 
odic  solutions  exhibiting  the  bifurcation  points  Ap/j  c;  -1.61  (analytical 
approximatiun:  Apo  —  —1.5)  and  Agp  —  —1,24  (analyticeil  approodmation: 
Aqp  ^  —1.35).  The  values  of  the  parameters  are  P  =  0.5,  b  =  S,  T  =  2000 
and  A  =  0.1. 

We  next  examine  a  case  of  positive  detuning  and  determine  the  bifurca¬ 
tion  diagram  in  terms  of  A  for  a  fixed  value  of  A  close  to  zero.  Figure  3  shows 
the  quasiperiodic  bifurcation  which  appears  soon  as  A  increases  from  zero. 
The  straight  lines  represent  the  maximum  and  the  minimum  of  the  small 
amplitude  multi- wave  mixing  regime  given  by  a  =  ±AA/(1  —  A^).  The  val¬ 
ues  of  the  parameters  are  P  =  0.5,6  =  5,  T  =  1000  and  A  =  0.25.  Note 
that  the  analytically  estimated  value  for  the  quasiperiodic  bifurcation  point, 
Agp  0.44,  compares  well  with  the  numerically  estimated  value,  Agp  ~ 

0.42. 
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